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ABSTRACT 


A  subsonic  jet  of  air  emanating  from  a  rectangular  nozzle  of  aspect 
ratio  16.7  has  been  investigated  using  hot-uire  anemometry.  The  purpose 
of  this  investigation  is  to  extend  our  knowledge  of  a  rectangular  jet 
from  incompressible  to  compressible  flow*  Four  exit  Mach  numbers*  rang¬ 
ing  from  0.18  to  0.8,  were  selected  for  the^studJ^T^  Results  of  the  hot¬ 
wire  measurements  show  that  a  compressible  subsonic  rectangular  jet 
behaves  similar  to  an  incompressible  jet.  It  exhibits  three  distinct 
regions  characterized  by  the  decay  of  the  mean  axial  velocity  along  the 
axis  of  the  jet.  These  three  regions  are  a  potential  core  region,  a 
two-dimensional  type  region,  and  an  axisymmetric  type  region.  The 
effect  of  increasing  the  Mach  number  of  the  rectangular  jet  was  found  to 
be  the  extension  of  the  two-dimensional  type  region,  and  the  jet  behaves 
more  nearly  like  a  two-dimensional  jet  at  higher  Mach  numbers.^  This 
effect  is  due  to  the  variation  in  the  growth  rate  of  the  jet  in  The 
plane  containing  the  long  dimension  of  the  nozzle  only.  The  Mach  number 
has  very  little  effect  on  the  spreading  in  the  plane  containing  the 
short  dimension  of  the  nozzle.  From  the  distributions  of  the  r.m.s. 
intensities  in  the  two  central  planes,  it  is  suggested  that  the  turbu¬ 
lence  structure  of  a  compressible  rectangular  jet  can  be  inferred  from 
that  of  an  incompressible  jet. 
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Chapter  1 


INTRODUCTION 

Rectangular  jets  have  been  used  in  a  wide  variety  of  engineering 
applications  such  as  the  augmentor  uing  ejectors  for  short-takeof f-and- 
landing  (STOL)  aircraft.  The  present  experiment  is  part  of  a  continuing 
effort  to  investigate  the  mixing  of  multiple  rectangular  jets  in 
thrust-augmenting  ejectors.  Previous  work  on  this  subject  has  been 
reported  by  Krothapalli  et  al.lll*  where  emphasis  was  on  the  incompres¬ 
sible  rectangular  jet.  The  present  work  concentrates  on  the  structure 
of  a  compressible  rectangular  jet  at  various  subsonic  Mach  numbers. 

The  previous  investigation  has  revealed  that  the  flow  field  of  an 
incompressible  rectangular  jet  can  be  characterized  by  the  presence  of 
three  distinct  regions  as  defined  by  the  decay  of  the  mean  axial  veloc¬ 
ity  along  the  centerline  of  the  jet.  Uc-  These  regions  may  be  classi¬ 
fied  as  follows:  a  potential  core  region,  where  Uc  is  constant;  a  two- 
dimensional  type  region,  where  Uc  decays  close  to  a  self-preserved 
two-dimensional  jet;  and  an  axisymmetric  type  region,  where  Uc  decays  at 
the  same  rate  as  a  self-preserved  axisymmetric  jet.  It  is  known  that 
the  extent  of  each  of  these  regions  is  a  function  of  the  nozzle  aspect 
ratio  and  inlet  geometry.  Detailed  mean  velocity  and  temperature  pro¬ 
files  as  well  as  some  turbulence  quantities  of  rectangular  jets  having 
different  aspect  ratios  and  nozzle  geometries  were  measured  using 


*  Numbers  in  brackets  designate  References  cited. 
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hot-uire  anemometry  by  Sfeir  [2,3*.  Recently.  Krothapalli,  Baganoff  and 
Karamcheti  [4]  reported  further  detailed  measurements  of  an 


incompressible  rectangular  jet.  Mhile  supporting  the  previous 
observations  (2>3f5-8],  these  results  furnished  more  turbulence  data. 

In  the  cases  of  thrust  augmenting  ejectors  and  jet  engines,  the  jet 
exit  velocities  are  so  high  that  the  flou  can  no  longer  be  considered  as 
being  incompressible.  However,  none  of  the  fore-mentioned  experiments 
were  conducted  at  velocities  high  enough  to  be  considered  compressible. 
It  is  the  purpose  of  this  investigation  to  study  the  flou  structure  of  a 
compressible  rectangular  jet  and  identify,  if  any.  the  effects  of  Mach 
number  on  the  jet  structure. 

In  light  of  past  experimental  observations  (1-8)  the  characteristics 
of  a  rectangular  jet  are  expected  to  depend  upon  the  following  parame¬ 
ters: 

•  aspect  ratio  of  the  nozzle 

•  inlet  geometry  of  the  nozzle 

•  Mach  number  at  the  nozzle  exit 

•  Reynolds  number  at  the  nozzle  exit 

•  type  of  exit  velocity  profile 

•  type  of  boundary  layer  at  the  nozzle  exit 

•  magnitude  of  the  turbulence  intensity  at  the  exit  plane  of  the 


-  2  - 


nozzle 


•  condition  of  the  ambient  medium  into  uhieh  the  jet  is  issuing 

In  the  present  experiment,  four  exit  Mach  numbers  (0.18,  0.3,  0.5,  and 
0.8)  were  chosen.  The  corresponding  Reynolds  numbers  based  on  the  uidth 
of  the  nozzle  extended  from  12,000  to  52,000.  A  nozzle  of  aspect  ratio 
16.7  (50  mm  long,  3  mm  wide)  was  used  in  this  experiment.  The  inlet 


geometry  of  the  nozzle  was  designed  to  obtain  a  low  turbulence  level  at 


Chapter  2 


APPARATUS,  INSTRUMENTATION,  AND  PROCEDURES 

2.1  AIR  SUPPLY  SYSTEM 

A  schematic  of  the  air  supply  system  is  shoun  in  figure  1.  A  Uorth- 
ington  Type  280  compressor  uas  used  in  conjunction  with  a  104  ft’  stor¬ 
age  tank  to  provide  a  blou-doun  type  air  supply.  The  compressor  is  able 
to  raise  the  pressure  in  the  storage  tank  at  a  rate  of  150  psi  per  hour, 
to  a  peak  pressure  of  2800  psi.  For  the  highest  mass  flou  rate  used  for 
this  investigation  (125  kg/hr),  a  fully  charged  tank  is  able  to  supply 
the  air  flou  for  a  period  of  about  tuo  hours. 

The  high  pressure  air  from  the  storage  tank  passes  through  a  shut-off 
valve  and  a  0.5  pm  filter  before  entering  a  tuo-stage  feedback  control 
regulator.  The  tuo-stage  regulator  is  able  to  drop  the  pressure  to  70 
psig  or  louer  and  maintain  a  constant  pressure  in  the  settling  chamber 
to  an  accuracy  better  than  ±0.1  cmHg,  which  corresponds  to  a  velocity 
variation  of  less  than  2  m/sec  at  the  nozzle  exit. 

Because  the  air  undergoes  a  throttling  process  at  the  regulator,  the 
flou  temperature  drops  considerably.  For  conditions  of  this  experiment, 
the  maximum  temperature  drop  is  about  10’C.  Provisions  uere  made  to 
raise  the  flou  total  temperature  to  that  of  ambient  so  that  the  total 
temperature  of  the  jet  remains  the  same  upon  mixing  uith  ambient  air. 
This  uas  done  by  a  heating  stage  which  consisted  of  sixteen  strip 
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heaters  (85  cm  long  and  A  cm  uide  each)  and  tuo  square  pipes  (2x2  in^ 
cross  section,  each  5  ft  long).  The  strip  heaters  were  mounted  on  the 
outside  walls  of  the  square  pipes  through  which  the  air  passes.  A 
picture  of  this  setup  is  shown  in  figure  2a.  Only  eight  of  the  sixteen 
heaters  were  used  in  this  experiment.  Each  heater,  which  can  produce 
7.5  kM  of  power,  was  connected  to  an  independent  switch.  A  schematic  of 
the  arrangement  is  shown  in  figure  2b.  The  heaters  are  controlled 
manually  by  turning  on  a  certain  number  of  heaters  while  monitoring  the 
temperature  in  the  settling  chamber  using  a  digital  thermometer  (Doric 
410A  Digital  Trendicator  with  iron-constantine  thermocouple).  For  the 
highest  flow  rate  used  in  this  experiment  (0.04  kg/sec),  the  total 
temperature  can  be  raised  up  to  40*^0  above  ambient  temperature.  The 
temperature  in  the  settling  chamber  was  maintained  at  a  constant 
temperature,  usually  at  room  temperature,  to  an  accuracy  of  about  tO.S^C 
over  the  duration  of  each  test. 

After  passing  through  the  heating  stage,  the  air  is  injected  into  a 
cylindrical  settling  chamber  through  a  spray  bar  (see  figure  1)  which 
helps  in  distributing  the  flow  uniformly.  The  settling  chamber  is 
1.75  m  in  length  and  0.6  m  in  diameter.  The  velocity  of  the  flow  in  the 
settling  chamber  is  about  0.14  m/sec  for  the  highest  flow  rate  tested. 
At  the  other  end  of  the  settling  chamber  the  air  is  passed  through  a 
bell  mouth  and  six  fine  mesh  screens  placed  5  cm  apart  in  the  exit  sec¬ 
tion  of  the  settling  chamber.  Following  the  screen  pack,  the  flow  pas¬ 
ses  through  an  adapter  which  connects  the  model  to  the  settling  chamber. 
The  adapter  was  made  to  provide  a  smooth  transition  from  a  circular 
cross  section  at  the  exit  of  the  settling  chamber  to  the  rectangular 


cross  section  of  the  model. 


This  experiment  made  use  of  the  same  air 


supply  system  as  the  one  used  by  Krothapalli  et  al.MI  except  for  the 
heating  stage.  More  detailed  information  on  the  facility  is  found  in 
that  report. 

2.2  NOZZLE  GEOMETRY  AND  COORDINATE  SYSTEM 

Figure  3a  shous  the  geometry  of  the  rectangular  nozzle  together  uith 
the  coordinate  system  used  in  this  study.  The  nozzle  is  50.0  mm  in 
length  (L)  and  3.0  mm  in  width  (D),  and  is  preceded  by  a  40  mm  long 
channel  having  the  same  cross  section.  The  nozzle  lip  has  a  thickness 
of  1.5  mm.  The  convergent  entrance  region  was  designed  to  provide 
smooth  contraction  uith  a  contraction  ratio  of  about  36.  The  nozzle 
used  is  the  center  lobe  of  a  multiple  nozzle  configuration*  as  shown  in 
figure  3b.  which  was  designed  to  study  the  mixing  of  multiple  rectangu¬ 
lar  jets  (1].  The  remaining  four  lobes  were  sealed  at  the  inlet  sec¬ 
tion.'  A  70-degree  triangular  shaped  ventilation  entry  between  the  lobes 
was  employed  to  ensure  a  smooth  induced  flow.  For  the  conditions 
tested,  top-hat  mean  velocity  profiles  and  low  turbulence  intensities 
were  obtained  at  the  nozzle  exit. 

A  Cartesian  coordinate  system  was  used  uith  its  origin  located  at  the 
center  of  the  nozzle  and  uith  the  X  axis  oriented  along  the  centerline 
of  the  jet.  The  Y  and  Z  axes  are  parallel  to  the  short  and  long  dimen¬ 
sions  of  the  nozzle,  respectively. 
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2.3  FLOM  VISUALIZATION  SYSTEfl 

A  conventional  Schlieren  system,  as  shown  schematically  in  figure  4, 
was  adopted  for  flow  visualization.  The  light  source  employed  is  a 
stroboscopic  flash  unit  (US  Scientific  Instruments  Type  3015  Strobo¬ 
scope)  which  can  be  triggered  either  manually  to  produce  a  single  flash 
or  as  a  stroboscope.  The  duration  of  each  flash  is  a  function  of  the 
light  intensity  and  can  be  varied  from  1.3  to  7.0  psec  at  five  discrete 
values.  Each  of  the  Schlieren  photographs  shown  in  this  report  was 
taken  using  a  single  flash  of  1.5  psec  duration.  The  collimating  and 
focusing  mirrors  are  spherical  mirrors.  Their  focal  lengths  (f)  and 
diameters  (dia.)  are  given  in  the  figure.  The  optics  were  arranged 
such  that  the  size  of  the  image  was  the  same  as  that  of  the  object. 
Polaroid  Type  57  high  speed  instant  film  (ASA  3000)  was  used  to  obtain 
the  photographs. 

2.4  HOT-HIRE  ANEMOMETRY  Mil  TRAV.ERSINQ 

The  velocity  measurements  were  made  using  a  linearized  constant-tem- 
perauture  hot-wire  anemometer.  The  hot-wire  probe  was  positioned  by  an 
orthogonal  three-axis  traversing  mechanism.  A  schematic  of  the  anemome¬ 
ter  with  the  signal  processing  instruments  and  the  traversing  mechanism 
is  shown  in  figure  5. 

A  standard  constant-temperature  bridge  DISA  55ni0  in  conjunction  with 
a  DISA  55P11  single  normal  wire  was  used  to  make  the  measurements.  The 
wire  is  a  platinum-coated  tungsten  wire  having  a  5  pm  diameter  and  a 
1  mm  length.  The  frequency  response  obtained  with  a  square  wave  test 
was  approximately  40  kHz.  All  hot-wire  measurements  were  made  with  the 
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Hire  oriented  parallel  to  the  Z  axis  and  perpendicular  to  the  X  and  Y 
axes.  With  the  objective  of  the  experiment  in  mind»  which  was  to 
explore  and  study  the  jet  behavior  at  high  velocities  in  terms  of  global 
properties  rather  than  detailed  information,  a  single  wire  was  used 
rather  than  an  X-uire  probe. 

To  linearize  the  hot-wire  signals,  a  DISA  551125  Linearizer  was  used. 
A  linearization  range  between  20  and  280  m/sec  was  achieved  in  this 
experiment.  Details  of  the  calibration  procedure  are  discussed  in  the 
next  section. 

The  signals  from  the  linearizer  were  sent  to  a  TSl  (Thermo-Systems 
Inc.)  Type  1076  Digital  Voltmeter  and  a  DISA  55D35  ft.m.s.  Voltmeter  for 
mean  and  r.m.s.  values.  The  integration  time  constant  can  be  set  in 
discrete  steps  from  0.1  to  100  seconds  on  both  voltmeters. 

A  traversing  mechanism  supplied  by  Velmex  Corp.  was  employed  to  posi¬ 
tion  the  hot-wire  probe.  It  has  three  translational  degrees  of  freedom. 
Each  of  the  three  axes  was  oriented  parallel  to  an  individual  coordinate 
of  the  Cartesian  system  chosen  (see  figure  3a).  The  vertical  axis  of 
the  mechanism  (parallel  to  the  X  axis)  was  equipped  with  a  Bodine  induc¬ 
tion  synchronous  motor  which  allowed  quick  positioning  of  the  probe  to  a 
required  downstream  location.  The  two  horizontal  axes  of  the  traversing 
mechanism  (parallel  to  the  Y  and  Z  axes)  were  guided  by  lead  screws  and 
driven  separately  by  two  stepper  motors.  The  stepper  motor  rotates  the 
lead  screw  through  a  fixed  angle  for  each  pulse  sent  by  the  control  cir¬ 
cuitry.  The  total  angular  rotations  or  the  total  distance  moved  can  be 
determined  by  keeping  track  of  the  number  of  pulses  sent  to  the  motor. 


The  minimum  positioning  increment  is  0.01  mm.  The  control  circuitry  can 
send  pulses  at  a  fixed  rate  as  selected  by  the  operator.  The  circuitry 
also  provides  an  analog  voltage  proportional  to  the  relative  position  of 
the  probe.  Uith  this  capability,  mean  and  r.m.s.  velocity  profiles  can 
be  recorded  in  a  continuous  mode.  This  uas  done  on  tuo  X-Y  recorders 
(HP  7004B  and  7046A).  The  X  axes  of  the  recorders  uere  connected  to  the 
analog  output  of  the  traverse  control  circuitry  while  the  Y  axes  uere 
connected  separately  to  the  outputs  of  TSI  and  DISA  voltmeters.  Having 
tested  different  voltmeter  time  constants  and  traversing  speeds,  it  uas 
found  that  the  combination  of  a  1  second  time  constant  and  a  traversing 
speed  of  approximately  10  mm/min  uas  adequate,  and  hence  used  in  this 
experiment. 

2.5  CALIBRATION  EQUIPMENT  AND  PROCEDURES 

The  output  voltage  of  a  constant-temperature  anemometer  (eta)  may  be 
expressed  in  the  following  form  [9]: 

Ecta^  =  A  +  B  (pU)"  (1) 

where  Ecfa  is  the  anemometer  output  voltage,  p  and  U  are  the  density  and 
velocity  of  the  fluid,  respectively,  at  the  point  where  the  hot-wire 
sensor  is  located,  and  A.  B  and  n  are  constants  to  be  determined  experi¬ 
mentally.  The  value  of  the  exponent  n  is  always  close  to  0.5.  Equation 
(1)  shows  that  the  anemometer  output  voltage  is  a  nonlinear  function  of 
the  mass  flow  rate  pD.  The  DISA  55n25  Linearizer  accomplishes  the  line¬ 
arization  by  means  of  an  analog  computer  having  a  transfer  function  of 
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the  following  form: 


E 1 i n  ”  Kl  (  Ecta 


A  )" 


(2) 


2 

uhere  Eiin>  Ki  and  m  are  the  tinearizer  output  voltage,  a  constant  and 
the  exponent,  respectively.  In  the  particular  case  of  the  DISA  55n25 
Linearizer  the  exponent  m  can  be  made  velocity  dependent.  m=m(0).  Sub¬ 
stitution  of  Eq.(l)  in  (2)  gives 

Eiin  =  Ki  B"  (pU)""  (3) 

In  the  case  of  incompressible  flows.  the  density  is  constant  and 
letting  m=1/n  a  direct  proportionality  between  Eiin  and  U  is  given 

Enn  =  K  U  (4) 

where  K  is  chosen  by  the  linearizer  gain  control. 

In  compressible  flow  the  density  p  can  be  expressed  as  a  function  of 
the  velocity  U  and  then  with  a  properly  chosen  m(U)  a  direct  relation¬ 
ship  between  the  linearizer  voltage  and  the  flow  velocity  can  also  be 
achieved.  A  cold  jet  exhausting  into  a  medium  at  rest  is  an  example  of 
such  a  flow  and  has  been  used  to  calibrate  the  hot-wire.  This  flow  was 
supplied  by  a  DISA  55090  Calibration  Equipment.  The  equipment  requires 
a  compressed  air  supply  of  200  psig  to  produce  a  free  air  jet.  at  low 
turbulence  level,  for  a  velocity  range  between  0.5  and  3Ui  m/sec.  The 
jet  velocity  with  correction  for  compressibility  is  available  as  an  ana¬ 
log  voltage  that  can  be  used  to  plot  the  calibration  curve  on  an  X-Y 
recorder. 
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Although  this  equipment  permits  rapid  and  efficient  calibration!  it 
has  no  control  of  the  flou  temperature  which  is  dependent  upon  the  air 
supplied.  Knowing  the  difference  between  the  calibration  temperature 
and  that  of  the  flow  being  measured,  one  could  use  analytical  formulas 
to  make  corrections  on  the  hot-wire  measurements.  However#  the  present 
experiment  used  a  different  approach  to  eliminate  the  error.  The  proce¬ 
dure  is  described  as  follows. 

After  linearization  is  achieved  on  the  DISA  calibrator#  the  hot-wire 
is  placed  at  the  center  of  the  exit  of  the  rectangular  nozzle.  Since 
the  flow  at  the  nozzle  exit  is  isentropic  and  uniform#  the  velocity  can 
be  calculated  from  the  total  pressure  and  temperature  by  the  isentropic 
flow  relations.  Keeping  the  total  temperature  the  same  as  that  of  the 
ambient  air  and  varying  the  total  pressure#  the  linearizer  gain  is 
adjusted  such  that  the  calibration  constant  K  equals  the  selected  value. 
This  procedure  was  repeated  before  each  test  run  and  the  temperature  was 
kept  the  same  during  the  test.  Figure  6  snows  the  resulting  linearized 
relation  between  the  mean  hot-wire  output  voltage  and  the  mean  flow 
velocity  for  the  range  between  20  and  280  m/sec  with  the  calibration 
constant  K  equal  to  1/28.  The  deviation  of  the  output  voltage  from  a 
straight  line  is  less  than  il%  for  velocities  higher  than  40  m/sec.  The 
eta  bridge  voltage  is  also  included  in  the  figure  to  show  the  nonlinear 
relation  between  the  mean  eta  voltage  E^ta  and  the  mean  velocity  U. 
Using  the  same  procedure#  different  velocity  ranges  were  chosen  for 
measurements  in  jets  at  different  exit  Mach  numbers.  For  each  calibra¬ 
tion  range#  the  maximum  output  voltage  of  the  linearizer  was  adjusted  to 
be  equal  to  10.0  volts  (maximum  linearizer  output)  for  the  best  resolu¬ 
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tion. 


2.6 


IkQM 


INTERPRETATION  (DF  HOT-»IRE  SIGNALS  m  COMPRESSIBLE 
For  measurements  in  a  "steady"  turbulent  jet,  the  linearizer  output 
voltage  Eiin  can  be  decomposed  into  a  time-averaged  mean  voltage  E  and  a 
fluctuation  voltage  e. 

E 1  i  n  -  E  +  e 


uhere  the  time  average  of  e  is  zero  by  definition.  E  is  directly  pro¬ 
portional  to  the  mean  velocity  U  as  the  calibration  shown  in  figure  6. 
The  fluctuation  voltage  e,  however,  is  dependent  upon  fluctuations  in 
mass  flow  rate,(pu)',  and  in  total  temperature.  To'.  Since  the  total 
temperatures  between  the  jet  flow  and  ambient  are  equal  in  the  present 
study,  it  may  be  assumed  that  the  total  temperature  fluctuation  To'  is 
negligible  and  e  becomes  a  measure  of  (pu)'  only. 

To  separate  u  from  (pu)',  Horstman  and  Rose  (101  suggested  an  analyt¬ 
ical  formula  with  an  assumption  that  <p'>  (r.m.s.  of  the  pressure  fluc¬ 
tuation)  is  negligible.  By  letting  <To'>  and  <p'>  be  zero,  the  formula 
given  by  Horstman  and  Rose  (10)  can  be  simplified  to 


<u>  1  <(pu)'> 


U  1+(y-1)M  pU 


(5) 


uhere  M  is  the  Mach  number  of  the  flow  at  the  sensor  location  and  y  the 
specific  heat  ratio.  By  assuming  the  flow  to  be  locally  isentropie,  the 
Mach  number  M  can  be  calculated  from  mean  velocity  U.  Using  the  value 
of  <e>/E  for  <(pu)'>/pU  in  equation  (5),  <u>/U  can  be  found.  This 
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method  uas  used  to  find  <u>/U  (denoted  as  u)  from  measured  <e>/E  (=  S) 
along  the  centerline  of  the  jet  for  exit  Mach  number  of  0.8.  A  compari¬ 
son  of  the  values  between  u  and  $  is  displayed  in  figure  7.  At  X/D=150> 
the  flow  velocity  is  about  50  m/sec  which  is  essentially  an  incompressi¬ 
ble  flow  but  the  formula  still  gives  about  5  per  cent  correction  for  u. 
Because  of  this  doubtful  result*  no  attempt  is  made  here  to  incorporate 
any  compressibility  corrections  to  the  velocity  fluctuations.  The  fluc¬ 
tuation  quantities  will  be  expressed  as  r.m.s.  of  fluctuation  voltage* 
<e>*  normalized  with  centerline  mean  voltage  Ee  at  the  station  in  ques¬ 
tion  and  expressed  as  e  which*  in  an  incompressible  flow*  is  equivalent 
to  u. 

2.7  TRANSONIC  EFFECT 

According  to  Liepmann  and  Roshko  Mil*  a  flow  is  called  transonic  if 
both  subsonic  and  supersonic  regions  are  present  in  the  flow  field.  By 
this  definition  a  jet  of  exit  nach  number  of  0.8  is  a  bona  fide  subsonic 
flow.  However*  as  soon  as  a  probe  is  inserted  in  the  flow*  it  may 
become  transonic  and  the  flow  characteristics  may  be  changed. 

To  ensure  that  such  an  effect  is  not  found  in  this  experiment*  the 
hot-wire  with  its  support  uas  placed  in  the  potential  core  of  the  rec¬ 
tangular  jet  at  the  highest  exit  Mach  number  tested  (0.8)  and  investi¬ 
gated  using  Schlieren  system.  Figures  8a  and  b  show  the  flow  with  and 
without  the  probe  in  place.  On  comparing  these  two  photographs*  no 
obvious  effect  of  the  probe  is  noticed. 
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Chapter  3 


RESULTS  AND  DISCUSSION 

In  this  chapter,  most  of  the  measurements  presented  are  for  exit  Mach 
numbers  (Me)  of  0.3>  0.5  and  0.8.  For  comparison  purposes>  uherever 
appropriate,  results  of  the  incompressible  jet  exiting  at  Me  equal  to 
0.18>  taken  from  a  previous  investigation  [4],  are  included. 

3. 1  INITIAL  CONDITIONS 

It  is  known  that  the  development  of  a  rectangular  jet  depends  on  the 
conditions  at  the  nozzle  exit.  In  the  present  investigation,  the  condi¬ 
tions  at  the  nozzle  exit  have  been  studied  using  the  hot-wire  and  pres¬ 
ented  in  figure  9.  The  mean  and  r.m.s.  hot-wire  voltages  along  the  Y 
and  Z  axes  at  the  nozzle  exit  are  recorded  for  Me  equal  to  0.3.  0.5  and 
0.8.  The  profiles  along  the  Y  axis  have  been  stretched  ten  times  as 
wide  as  those  along  the  Z  axis  in  order  to  show  the  details  in  the  short 
dimension  of  the  jet.  For  optimum  resolution  of  the  profile,  different 
calibration  constants  K  (defined  as  E=KU)  were  chosen  (by  varying  the 
linearizer  gain)  to  make  the  measurements  in  the  jet  for  different  Mach 
numbers. 

Along  the  Y  axis,  top-hat  mean  voltage  profiles  were  obtained  for  all 
three  exit  Mach  numbers.  The  r.m.s.  profiles  along  the  Y  axis  show 
peaks  in  the  shear  layers  on  either  side  of  the  jet.  The  r.m.s.  inten¬ 
sity  (<e>/Ec)  in  the  shear  layer  increases  with  Mach  number.  At  the 
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center  of  the  nozzle  exiti  the  r.m.s.  intensities  are  0.3%,  0.3%  and 


I 


0.45%  for  M«  equal  to  0.3,  0.5  and  0.8,  respectively. 

Along  the  Z  axis,  top-hat  mean  voltage  profiles  uere  also  obtained 
except  for  slight  variations  close  to  the  jet  boundaries  (approximately 
10  mm  inside  the  boundary).  The  r.m.s.  profiles  shou  spikes  at  the  jet 
boundaries  and  the  r.m.s.  intensity  (<e>/Ec)  of  the  peak  increases  with 
Me.  At  the  Z  coordinate  corresponding  to  the  location  uhere  the  mean 
voltage  profile  shows  non-uniformity,  the  r.m.s.  profile  has  a  local 
maximum  value  which  decreases  in  magnitude  as  Me  increases.  It  is  sug¬ 
gested  that  this  local  turbulence  may  be  a  result  of  the  secondary  flow 
generated  in  a  rectangular  channel  ahead  of  the  nozzle  exit  (see 
Schlichting  [121). 

3.2  SCHLIEREN  PHOTOGRAPHS 

To  study  the  jet  qualitatively,  Schlieren  photographs  of  the  jet  with 
various  exit  Mach  numbers  uere  taken.  Figures  10a, b,c  shou  pictures  of 
the  jet  in  the  plane  containing  the  short  dimension  (X,Y  plane)  of  the 
nozzle  for  equal  to  0.3,  0.5  and  0.8.  These  pictures  shou  the  flou 
field  of  the  jet  from  the  nozzle  exit  to  a  downstream  distance  of  about 
30  nozzle  widths.  No  significant  variation  in  the  spreading  rate  of  the 
jet  with  Mach  number  is  noticed  in  these  photographs. 

Figure  11  shows  the  jet  for  M«  equal  to  0.8  in  the  plane  containing 
the  long  dimension  (X,Z  plane)  of  the  nozzle.  Mithin  the  field  of 
observation,  very  little  spreading  of  the  jet  is  observed,  indicating 
that  the  flou  is  nearly  tuo-d imensiona) .  Photographs  of  the  jet  for 
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equal  to  0.3  and  0.5  are  not  included  because  of  their  poor  contrast; 
they  show  similar  features  as  that  for  Me  equal  to  0.8. 


3.3  MEAN  VELOCITY  FIELD 

The  variation  of  the  (squared)  normalized  mean  velocity  along  the 
centerline  of  the  jet  for  four  exit  Mach  numbers  is  shoun  in  figure  12. 
The  centerline  mean  velocity  is  normalized  with  respect  to  the  mean 
velocity  at  the  center  of  the  nozzle  exit.  The  three  regions  of  an 
incompressible  rectangular  jet  as  described  in  Chapter  1  can  be  identi¬ 
fied  in  this  figure.  The  first  region  is  the  potential  core  uhere  the 
centerline  mean  velocity  Uc  is  constant.  It  extends  from  the  nozzle 
exit  to  a  dounstream  location  of  about  5D.  In  the  second  region,  so 
called  tuo-dimensional  type  region,  the  decay  rate  is  close  to  that  of  a 
tuo-dimensional  jet,  (Uc^Uo)^*  X*'.  In  the  third  region,  the  jet  cen¬ 
terline  mean  velocity  decays  similar  to  an  axisymmetr-'r  ;et, 
(Uc/Uo)*a  X**,  hence  called  axisymmetric  type  region.  It  is  Jtserved 
that  the  transition  from  the  second  to  the  third  region  seems  to  occur 
slightly  farther  dounstream  as  M,  increases.  It  appears  that  the  jet 
decay  rate  has  a  fairly  ueak  dependence  on  the  Mach  number. 

Figures  13a, b,c  shou  the  distribution  of  mean  velocity  (U)  in  the 
central  X,Y  plane  at  different  dounstream  locations,  ranging  from  20  to 
100  nozzle  uidths(D)  for  various  exit  Mach  numbers.  The  velocity  U  is 
normalized  with  respect  to  Uc  at  the  corresponding  dounstream  location 
uhile  the  distance  Y  is  normalized  by  the  distance  X  to  the  station  in 
question.  The  similarity  profile  for  a  plane  turbulent  jet,  from 
Heskestad  113],  is  shoun  for  comparison  in  each  of  the  plots.  Uithin 
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the  limits  of  error  for  the  experiment*  at  each  Mach  number  the  profiles 
are  geometrically  similar.  The  shape  of  the  similarity  profile  is 
nearly  independent  of  the  exit  Mach  number  and  similar  to  that  of  a 
two-dimensional  jet.  The  similarity  profile  for  an  incompressible 
rectangular  jet  [4]  is  included  in  figure  13c.  It  appears  that  in  the 
central  X>Y  plane,  compressibility  does  not  introduce  a  siginificant 
effect  on  the  self-similarity  of  the  mean  velocity  profile. 

Figures  14a, b,c  show  normalized  mean  velocity  profiles  in  the  central 
X,Z  plane  at  five  downstream  locations  for  various  exit  Mach  numbers. 
At  each  corresponding  downstream  location,  the  profile  shows  only  slight 
variation  with  Me.  At  X=2QD,  a  saddle  shape  profile  is  observed  in  the 
jet  of  Me  equal  to  0.3.  This  type  of  profile  was  also  observed  in  the 
two-dimensional  type  region  of  incompressible  rectangular  jets  11-71. 
As  the  exit  Mach  number  increases  to  0.5  and  0.8,  the  saddle  shape 
becomes  less  pronounced.  These  mean  velocity  profiles  suggest  that  in 
the  two-dimensional  type  region  the  jet  behaves  more  nearly  like  a  two- 
dimensional  jet  at  higher  subsonic  Mach  numbers.  The  similarity  profile 
in  the  central  X,Y  plane  (figure  13)  is  also  included  in  figure  14.  On 
comparing  the  profiles  between  the  two  central  planes,  it  is  observed 
that  the  location  where  the  profiles  in  both  planes  are  identical  seems 
to  move  slightly  downstream  with  increasing  Mach  number.  As  will  be 
seen  later,  this  location  corresponds  to  the  "crossover  point".  Down¬ 
stream  from  this  point,  the  shape  of  the  profile  approaches  that  of  an 
axisymmetric  jet  as  observed  in  the  incompressible  jet  (41.  The  profile 
of  a  self-preserved  axisymmetric  jet  measured  by  Uygnanski  and  Fiedler 
1141  is  included  in  figure  14a.  As  compared  to  an  incompressible  jet, 
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the  variations  noted  above  are  not  significant  enough  to  have  a  major 
influence  on  the  overall  development  of  the  jet. 


The  growth  of  the  jet  with  downstream  distance  in  the  central  XfY  and 
XpZ  planes  can  be  represented  by  the  half-velocity  widths  Y1/2  and  Z^/^, 
respectively,  as  shown  in  figure  15  for  different  exit  Mach  numbers. 
The  half-velocity  width  is  defined  as  the  distance  from  the  centerline 
of  the  jet  to  the  point  where  the  axial  mean  velocity  is  equal  to  one- 
half  of  its  centerline  value.  The  jet  in  the  X.Y  plane  spreads  linearly 
with  X  and  the  locus  of  the  half-velocity  points  is  given  by 


Yi/2  =  k  (  X  -  Xo) 

where  k  is  the  spreading  rate  and  X©  the  virtual  origin.  The  values  of 
k  and  X©  are  unchanged  for  M©  varying  from  0.3  to  O.S  (k=0.11>  X©=-2.00) 
and  show  very  little  difference  from  those  of  an  incompressible  jet  at 
M,  equal  to  0.18  (k=0.109.  X©='2.5D)  141.  The  exit  Mach  number  seems  to 
have  very  little  effect  on  the  growth  of  the  jet  in  the  X.Y  plane. 

However,  in  the  central  X.Z  plane  the  variation  of  the  half-velocity 
width  with  X  shows  some  changes  with  exit  Mach  number.  Tor  M©  equal  to 
0.3.  0.5  and  0.8.  Zi/2  varies  only  slightly  with  X.  But  in  the  incom¬ 
pressible  jet  (M©=0.18)  [41.  Z1/2  first  decreases  with  downstream  loca¬ 
tion  X  and  then  increases.  At  some  intermediate  location  the  half-ve¬ 
locity  widths  in  the  two  central  planes  cross  over.  This  location 
corresponds  to  the  onset  of  the  axisymmetric  type  region.  The  distance 
from  the  nozzle  exit  to  the  crossover  point  along  the  X  axis  is  denoted 
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as  Xe<  It  is  noted  that  Xc  increases  with  increasing  Mach  number,  uhich 
is  solely  due  to  the  variation  of  Zi/j  uith  the  exit  Mach  number. 

3.4  S.ri.S.  INTENSITY 

The  r.m.s.  values  of  hot-wire  voltage  fluctuations  on  the  centerline 
of  the  jet  normalized  uith  the  local  mean,  e  (=<e>/Ee)«  for  four  differ¬ 
ent  exit  Mach  numbers  are  shown  in  figure  16.  The  variations  of  e  along 
the  jet  centerline  for  all  four  Mach  numbers  follow  the  same  pattern. 
The  magnitude  of  e  increases  sharply  in  the  potential  core  region.  In 
the  two  dimensional  type  region,  e  either  decreases  or  maintains  a  con¬ 
stant  value  uith  X  and  then  increases  gradually.  The  gradual  increase 
in  e  is  continued  in  the  axisymmetric  type  region.  For  X  greater  than 
about  600.  the  magnitude  of  S  at  a  given  downstream  location  seems  to 
decrease  uith  increasing  M*. 

The  behavior  of  $  described  above  for  X  less  than  300  is  very  similar 
to  that  of  a  rectangular  jet  uith  laminar  boundary  layers  at  the  nozzle 
exit  as  reported  by  Hill.  Jenkins  and  Gilbert  [15]  (aspect  ratio  10.5). 
They  also  found  that.  when  the  boundary  layers  are  turbulent.  $ 
increases  monotonical ly  uith  X.  On  comparing  the  variations  of  9  in 
figure  16  uith  the  results  in  Ref. 15.  it  may  be  concluded  that  the 
boundary  layers  in  the  present  experiment  are  laminar  for  all  four  exit 
Mach  numbers  tested.  However,  the  turbulence  intensity  in  the  boundary 
layer  at  the  nozzle  exit  for  M*  equal  to  0.8  was  found  to  be  about  10)1 
(see  figure  9).  uhich  led  us  to  believe  that  the  boundary  layer  might  be 
turbulent.  Due  to  the  small  size  of  the  nozzle,  no  measurements  could 
be  made  in  the  boundary  layer  to  justify  either  one  of  the  above  state¬ 
ments. 
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Profiles  of  e  in  the  central  )(#Y  plane  at  different  dounstream  loca¬ 
tions  are  shoMn  in  figures  17a>b»c,  for  various  exit  Mach  numbers.  A 
distinct  saddle  shape  is  observed  in  the  profiles  in  both  tuo-dimen- 
sional  and  axisymmetric  type  regions.  The  e  profiles  indicate  geometri¬ 
cal  similarity  only  in  the  tHo-d'<fflensional  type  region  of  the  jet 
(200  <  X  <  800).  This  uos  also  observed  in  the  incompressible  jet  141. 
For  X  greater  than  800  (i.e.  in  the  so-called  axisymmetric  type  region), 
the  magnitude  of  S  near  the  jet  centerline  increases  and  the  profiles 
approach  that  of  an  axisymmetric  jet  (see  figure  A  of  [14]).  This 
behavior  is  believed  to  be  a  result  of  the  merging  of  the  tuo  shear  lay¬ 
ers  separated  by  the  long  dimension  of  the  nozzle. 

Profiles  of  e  in  the  central  X,Z  plane  are  shown  in  figures  18a, b,c 
for  different  M«.  On  comparing  these  profiles  with  the  corresponding 
mean  velocity  profiles  in  figure  14,  the  point  of  maximum  turbulence 
intensity  coincides  with  the  point  where  the  velocity  gradient  (bU/M) 
is  maximum.  The  development  of  the  profiles  with  downstream  distance  in 
the  present  case  is  quite  similar  to  those  found  in  an  incompressible 
jet  [4J. 

3.5  AN  OVERALL  VIEH  OF  THE  FJ.014  FIELD 

In  order  to  provide  an  overall  picture  of  the  development  of  a  rec¬ 
tangular  jet,  isometric  plots  were  made  from  mean  and  r.m.s.  hot-wire 
voltages  of  the  jet  with  I1«  equal  to  0.8  at  three  downstream  locations 
(ISO,  400  and  900),  as  shown  in  figure  19.  Due  to  the  symmetry  of  the 
jet  about  the  central  X,Y  plane,  only  one  half  of  the  profiles  are  dis¬ 
played.  At  X  equal  to  150,  the  jet  is  nearly  two  dimensional  with  a 
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typical  saddle  shape  r.m.s.  profile  in  the  X,Y  plane.  Moving  dounstream 
to  400>  the  jet  spreads  in  the  X>Y  plane  but  with  no  apparent  grouth  in 
the  X>Z  plane.  At  X  equal  to  900>  both  mean  and  r.m.s.  profiles  become 
almost  axisymmetric  about  the  X  axis  uith  a  ueak  saddle  shape 
r.m.s.  profile:  and  the  jet  flow  seems  to  have  lost  its  memory  of  the 
nozzle  geometry. 
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Chapter  4 


CONCLUSIONS 

Based  on  the  results  discussed  in  Chapter  3  the  following  conclusions 
are  drawn.  In  the  case  of  a  subsonic  rectangular  jet/  the  general  char¬ 
acteristics  of  the  flow  field  do  not  change  significantly  with  exit  Mach 
number.  Similar  to  the  case  of  an  incompressible  rectangular  jet#  the 
mean  flow  field  of  the  jet  is  characterized  by  the  presence  of  three 
distinct  regions  as  defined  by  the  decay  of  the  centerline  mean  veloc¬ 
ity.  These  three  regions  are:  a  potential  core  region,  a  two-dimen¬ 
sional  type  region,  and  an  axisymmetric  type  region.  The  onset  of  the 
second  region  where  the  shear  layers  separated  by  the  short  dimension  of 
the  nozzle  meet  appears  to  be  independent  of  the  jet  exit  Mach  number. 
However,  the  location  where  the  shear  layers  separated  by  the  long 
dimension  of  the  nozzle  meet  seems  to  occur  farther  downstream  with 
increasing  Mach  number. 

From  the  measurements  in  the  central  X.Y  plane  of  the  jet.  the  fol¬ 
lowing  conclusions  can  be  drawn.  The  mean  velocity  profiles  are  geome¬ 
trically  similar  with  the  shape  of  the  profile  almost  identical  to  that 
of  an  incompressible  two-dimensional  jet.  The  half-velocity  width  of 
the  jet  varies  linearly  with  downstream  distance,  with  its  slope  and 
virtual  origin  nearly  independent  of  the  exit  Mach  number.  The  r.m.s. 
profiles  show  self-similarity  only  in  the  two-dimensional  type  region 
with  their  shape  being  independent  of  the  Mach  number.  In  the 


ax i symmetric  type  region 


no  geometrical  similarity  uas  observed  in 


those  profiles. 

In  the  central  X,Z  plane,  the  saddle  shape  mean  velocity  profile 
observed  in  the  tuo-dimensional  type  region  of  an  incompressible  jet 
becomes  less  pronounced  as  the  exit  Mach  number  increases.  This  obser¬ 
vation  along  with  the  results  in  the  X.Y  plane  suggests  that  the  jet 
behaves  more  nearly  like  a  tuo-dimensional  jet  at  higher  Mach  numbers 
for  downstream  locations  less  than  100D. 

From  these  results,  it  appears  that  for  a  given  nozzle  geometry  the 
overall  properties  of  a  subsonic  compressible  rectangular  jet  are  quite 
similar  to  that  of  an  incompressible  jet.  It  may  be  suggested  that  the 
turbulence  structure  of  a  compressible  rectangular  jet  can  be  inferred 
from  the  measurements  of  an  incompressible  rectangular  jet. 
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Figure  5  Signal  processing  instrumentation  and  probe  positioning  setup 


Figure  6  Calibration  curves  for  the  hot-wire 


Figure  7  Measured  r.m.s.  intensity  8  and  corrected  turbulence 
intensity  Q  along  the  centerline  of  the  rectangular 


Figure  9  Mean  and  r.m.s.  profiles  of  hot-wire  voltage  along  Y  and  Z  axes 
at  the  nozzle  exit  for  M_=  0.3,  0.5  and  0.8 


different  exit  Mach  numbers  (M 


Figure  12  The  decay  of  normaliz 
for  four  exit  Mach  nui 


Figure  13  (continue) 
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Figure  I3  (continue) 


Figure  14  (continue) 
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Lgure  16  Variation  of  turbulence  intensity  along  the  centerline 
of  the  jet  for  different  exit  Mach  numbers 


Distribution  of  hot-wire  voltage  fluctuations  in  the  central  X,Z  plane 
for  three  exit  Mach  numbers 


